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CT 1 Abstract 
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We study the generic transport of slice energy between the scalar field generated 
by the conformal transformation of higher-order gravity theories and the matter 
component. We give precise relations for this exchange in the cases of dust and 
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perfect fluids. We show that, unless we are in a stationary spacetime where slice 
energy is always conserved, in non-stationary situations contributions to the total 
slice energy depend on whether or not test matter follows geodesies in both frame 
representations of the dynamics, that is on whether or not the two conformally 
related frames are physically indistinguishable. 



1 Introduction 



With recent advances in observational cosmology [T] f(R) theories (and the closely re- 
lated family of scalar-tensor ones) have in the last few years regained much attention 
both in cosmology [2] and in other contexts [3] . Perhaps the most economical and conve- 
nient way to study such modified gravity theories is through their well-known conformal 
relation to general relativity, a technique developed in the eighties by different groups 

In this conformal method the /(-R)-vacuum equations on a spacetime (V,g) (the so- 
called Jordan frame) are transformed via a conformal transformation of the type g = e^g, 
where is a function of df/dR, to become on the conformally related spacetime (V, g) 
(the so-called Einstein frame) Einstein equations for the metric g with the scalar field 
having a self-interacting potential that depends on the function f(R) and its first 
derivatives. This technique provides a refreshing way to view the /(-R)-vacuum theory 
as a unified theory of gravitation (described in the Einstein frame by the metric g) and 
the scalar field 0, that IS ClS db theory uniting general relativity and the (lagrangian) 
theory of the scalar field. In the Jordan frame only one single geometric object appears, 
the 'metric' g, and the conformal transformation then serves as a tool to 'fragment' g 
into its two pieces in the Einstein frame, namely the gravitational field g and the scalar 
field 0. 

In a higher-order gravity theory with matter we have the following different pieces of 
information involved in the conformal transformation: 

• gravity, the field g or the conformally related field g 

• the scalar field 0. 

• the various matter fields if) which couple non-minimally to g and to 0, or their 
conformal transform ip which couples minimally to g but is not coupled to 0. 

The conformal transformation then relates the different pieces of matter and spacetime 
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geometry and describes the interaction between the components listed above in the 
context of f(R) theories. 

We describe in this paper how interactions of this sort lead to an exchange of slice 
energy between the various fields and spacetime geometry. For more varied interactions 
and energy transfer models of interest to cosmological situations see jH [7] and references 
therein. 

The plan of this paper is as follows. The next Section is preliminary and includes 
two simple applications (Corollaries 2.1, 2.2) of the basic properties (Theorems 2.1, 2.2) 
of the slice energy. These basic properties and also their proofs are included here for 
easy reference and also to establish notation. Their applications deal with the simpler 
case of f(R) theories in vacuum. Section 3 is the heart of this paper. There, we find the 
general (slice) energy transport equation and study in detail its application to the case 
of a general f(R) theory coupled to matter, where by matter we mean a general perfect 
fluid-scalar field system. There are many properties analysed here but a particularly 
important one deals with conditions under which the choice of 'physical' metric influences 
the net contribution to slice energy of the system. We conclude with a discussion of these 
results in Section 4. 

2 Energy on a slice 

Consider a time-oriented spacetime (V, g) with V = M. x K, where M. is a smooth 
manifold of dimension n, g a spacetime metric and the spatial slices JA t (= M. x {t}) 
are spacelike submanifolds endowed with the time-dependent spatial metric g t . (In the 
following, Greek indices run from to n, while Latin indices run from 1 to n. We also 

assume that the metric signature is (H — • ).) On (V,g) we consider a family of 

matterfields denoted collectively as if), assume that the field if) arises from a lagrangian 
density which we denote by L and denote the stress tensor of the field if) by T{if>). 

For X any causal vectorfield of V we define the energy-momentum vector P of a stress 
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tensor T relative to X to be 

P 13 = X a T af3 . (2.1) 
The energy on the slice M. t with respect to X is defined by the integral (when it exists) 



Et= / P a n a dnt, (2.2) 

JMt 

where n is the unit normal to M. t and d/i t is the volume element with respect to the 
spatial metric g t . We call P a n a the energy density. Assume that X and T are smooth. 
Then we have 

V Q P Q = V a {X p T^) = V a XpT a/3 + X p V a T^ 

or, equivalently, 

V Q P Q = l -T^(y a X p + V p X a ) + X p V a T^. (2.3) 

Thus, if /C C V is a compact domain with smooth boundary <9/C, it follows from Stokes' 
theorem that 

/ V a P a d/2= [ P a n a da, (2.4) 

JK JdK 

where d\i is the volume element of V and da that of <9/C, and so we find 

/ P a n a da= l - [ T al3 (V a X /3 + V (3 X a )dfi+ [ X^ a T^dfi. (2.5) 

JdK 1 JK JK 

Hence, when M. is compact or the field falls off appropriately at infinity, on the spacetime 
slab V = £ x [t , £ C M, and with T having support on V we have the following 
relation for the energies on the two end-slices 



P a n a dfi tl - ! P a n a dfi t0 = If' f T a P(V a X p + V p X a )dLi 

JM t0 1 J t JMt 

+ f 1 [ XpVaT^dfi (2.6) 

Jtn JMt 



or 

rti 



E tl - E t0 = \ ! 1 [ T^{V a X p + V p X a )dfi + ( 1 / XfF a r#dn. (2.7) 

Z J t JM t J t JMt 



3 



We therefore see that when X is a Killing vectorfield the first term on the right-hand-side 
of Eq. (12.71) is zero and so we have 

E tl - E t0 = f 1 f X p V a T^dfi. (2.8) 

Jto J M t 

Thus we have shown the following result (cf. [SJ, p. 87-88). 

Theorem 2.1 When X is a Killing vectorfield and the field is conserved, i.e., V a T al3 = 
0, we have 

E h =E t0 . (2.9) 

This means that, when the energy- momentum tensor of a field is conserved, the same is 
true for its slice energy relative to a Killing vectorfield as a function of time. 

In the next Section we pay particular attention to the case for which the field is a 
matter field if) interacting with a scalar field <fi with potential V ((/)). We take the scalar 
field lagrangian density to be 

L = -^d a <pdp<P + V(<P). (2.10) 
Then the energy-momentum tensor of (f) is 

T^(0) = Sf^t - ^(«9 a 0(9a0 - 2V{<t>)), (2.11) 
and we have the following result. 

Theorem 2.2 The energy density P a n a of the scalar field <p with potential V(<f>) is 
positive when V(<f>) > 0. 

Proof. The proof, which we give here for a general scalar field potential V(4>), is a 
direct adaptation with slight modifications of that found in [8], p. 88, for a power-law 
potential. Using Eq. (12.111) we calculate 

P a n a = ~X a n a d X (pd X (p + X a d u ^nPd^ + X a n a V{$). (2.12) 
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We define the quadratic form 

7 Am = - g ^x a n a + (XV + X"n A ) (2.13) 

and then we find that 

P a n a = ~X a n a g^d X (f>d^ + |(XV + X»n x )d x <pd^ + X a n a V(<P). (2.14) 
This means that 

P a n a = ^d^d^ + X a n a V -((f)). (2.15) 

Since M t is a t = const, hypersurface, we can choose coordinates such that X° = 1, X 1 = 
0,ni = 0. Then n = (# 00 )~ 1/2 , n l = g i0 (g 00 )~ 1/2 and so it follows that the quadratic 
form 7 is positive definite, 

7 00 = 2 00 n , 7*° = i ij = -g^no, (2.16) 

(recall signature of g^ is ( • — )). We then find that 

P a n a = X - (7 OO 9o09 o + 2 7 iO ^09 o + l ij d4d j( t>) + n V((f>) 
= \ (n g 00 <p 2 - g^nodrfdrf) + n Q V{4>) 

and therefore we conclude that the energy density P a n a is positive whenever V(<f>) > 0. 
This concludes the proof. ■ 

To end this Section, for the following simple application of the preceding developments 
we restrict attention to n = 4 spacetime dimensions although everything we do becomes 
valid with minor modifications to arbitrary n. The following notation for conformally 
related quantities is used: Let g and g be two conformal metrics, g = tt 2 g, on the 
manifold V. This means that in two orthonormal moving frames, 6 a and 6 a , the two 
conformal metrics satisfy 

g = Vafi e a e fi , g = r ]a p6 a 6 /3 and 9 a — Q _1 8 a , (2.17) 
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with r] a/ 3 = diag(+, — •••—) being the flat metric. Setting Q 2 = we see that 9 a = 
t-^l 2 B a and obviously 9 a = e<^ /2 # a . The same rules are true for any 1-form or vectorfield 
on V. Consider now the / (R)-vacuum equations, 

L aP = fR aP - X -g afi f - VaVpf + g a p □ g f = 0, (2.18) 

where the left hand side satisfies the conservation identities (cf. [9], p. 140) 

V a L afB = 0. (2.19) 

Then we conformally transform from (V,g) to the Einstein frame (V,g), according to 
the prescription given in [1], that is, we set 

= lnf, (2.20) 

to obtain the Einstein equations with a scalar field 'matter source' of potential V{4>) = 
(l/2)(f')~ 2 (Rf — /) and energy-momentum tensor given by Eq. (12. IIP : 

G aP = f ttP {<l>). (2.21) 

In this case we conclude that the field (ft is conserved, i.e., 

V a f af3 ((j)) = 0, (2.22) 

and, since 

V a f a/3 ((f)) = #V(V a S°0 + V), (2.23) 
we find that the 0-field is a scalar field satisfying the wave equation 

V a d a (j) + V = 0. (2.24) 
Further from Theorem 12. II we have the following result. 

Corollary 2.1 The slice energy of the scalar field </> generated by the conformal trans- 
formation 112. 20\) to the Einstein frame of the /(i?) -vacuum equations 112. 18\) relative to 
a Killing vectorfield of g is conserved, i.e., 

E t {4>) = / P a h a dfi t = const, (2.25) 

with dflt being the volume element of g t . 
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Secondly from Theorem 12.21 we have: 

Corollary 2.2 For all f(R)-vacuum theories A2.18\) with a positive potential in the Ein- 
stein frame the energy density P a h a of <fi is positive. 

Examples of theories in the last Corollary include, for instance, the choice f(R) = 
R + aR 2 ,a > 0. 

3 /(i?)-matter systems 

Suppose now that we start by coupling a matter field tj) to the geometry in (V, g) via the 
/(-R)-matter field equations 



Because of the conservation identities (12.191) . the field ip satisfies the conservation laws 



Then, if we conformally transform from (V, g) to the Einstein frame (V, g) according to 
(I2.20p . in place of equations (I2.2ip we obtain 



f'RaP ~ ~9apf - VaV/3/' + g a p □ g f = T a p(lp). 



(3.1) 



V Q T Q/3 (V>) = 0. 



(3.2) 



G af3 = T a p((f>) + T a p(if)) 



(3.3) 



where now the whole tensor in the right-hand-side is conserved, namely 




(3.4) 



but the two components are not conserved separately, that is 



V Q T a/3 (0) + 



(3.5) 



and 



V a f a P$) + 0, 



(3.6) 
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unless the conservation equations (13.21) for the field ip are conformally invariant (condi- 
tions for this are given in [10], p. 448). This result (already given in Ref. [11]) indicates 
that in higher order gravity theories there must be a generic, nontrivial — ip interaction 
between the matter field ip and the 0-field, and an associated exchange of energy between 
4> and ip. 

Writing Eq. (12. 7p for the scalar field <fi and substituting for the last term in the right- 
hand-side from Eq. (13. 4p we find the general energy transport equation in the Einstein 
frame, 



E tl (<P)-E tQ (<p) = - / f aP ((f>)(V a X +V p X a )dfl- / I^VX^)^, (3.7) 



with dft being the volume element of g. 

This result is only symbolic and has to be augmented by precise equations satisfied 
by the fields. Since the stress tensor of the 0-field is not separately conserved, it follows 
that the 0-field will not satisfy the usual scalar wave equation, □g0 + y'(0) = 0, but this 
equation will in general contain new terms. Similarly for the 'ordinary matter' ?/>-field, 
whatever its form (scalar field, Maxwell, a fluid etc), its field equations have new terms 
indicating the — ip interaction and associated energy exchange. For instance, if the ip 
field is another scalar field, then the equations satisfied by its conformal transform, ip, 
have a general form of the type O0 + h(<f)) d a (f) d a ip = 0, where h(4>) is a smooth function 
of 4> (often exponential). 

To study this interaction and the associated energy exchange between and ip more 
closely we give some concrete examples. Suppose firstly that -0 is a dust cloud on (V, g) 
with 4-velocity V a and stress tensor 





T a p, dust — pV a Vj3 



(3.8) 



satisfying the f(R)-dust equations in the Jordan frame, namely 



f'RafS - -g a pf ~ V Q V^/' + g aP □ g f = P V a Vp 



(3.9) 



Then 



v a { P v a vP) = 0, 



(3.10) 



8 



and it is obvious that here the dust streamlines are geodesies, that is V a is the tangent 
vectorfield to the geodesies. After the conformal transformation we find 

& a p = f aP (<i>)+pV a Vp, (3.11) 

with T a p((j)) given by (12.111) with tildes where appropriate, and with 

V a = e-*PV a , p = e-^p. (3.12) 

(We have set p = f2~ 4 p and since Vt 2 = e^, f2~ 4 = e~ 2( ^.) 

What is the field equation satisfied by the scalar field 0? From Eq. (13. lip the 
divergence of the stress tensor of is minus that of the dust, but 

V a (pV a V f) ) = V a (pV a V?) + A^pVW? + A^pV a V\ (3.13) 

where 

A% = \ (Spd 7 <t> + S^dp<f> - g^ ( g aS ds<P) . (3.14) 

From these equations and Eq. (13.101) we deduce the modified scalar field equation in the 
form 

d p <p{U<p + V) + -pV a V^d a (j) - ]-pd p 4> = 0. (3.15) 
Another way to derive the scalar field equation is as follows. Since 

V a T d f st = V a {~ P V a VP) = V^V a {pV a )+p{V a V^)V a , (3.16) 

and, since VpV? = 1, if we multiply Eq. (T2T231 by V p and use the fact that the divergence 
of the right hand side of Eq. A3. 11 j) is zero to arrive at the following equation for the 
scalar field <fi in the Einstein frame, namely, 

+ V) + V? V a {p V a ) + ~pV a V a V p = 0. (3.17) 

Recalling that dust matter follows geodesies on the original Jordan frame, V a 'V a V^ = 0, 
and taking it as a working hypothesis that the same is true in the conformally related 
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Einstein frame, we find that the last two terms in this equation are equal to the last two 
terms in Eq. (13 . 1 5[) and so we conclude that Eq. f!3.17|) provides an equivalent form of 
Eq. <^m>- 

We note that only in the very special case where we impose the constraint 



Vp = d p </>, 



(3.18) 



which implies some sort of 'alignment' between the dust component and the scalar field, 
does the scalar field equation (13.151) becomes the standard one, namely, 



D6 + V = 0. 



(3.19) 



We now study the behaviour of the total slice energy of the system comprised of 
and the dust component. We choose V = X so that 



P a n a = X p n aP V a V p = pV a n a . 



(3.20) 



Hence, applying Stokes' theorem we obtain 



K JdK. 



Therefore Eq. (13. 7\\ becomes 

! rti 



E t M - E ta {4>) 



f aP {4>){v a % + VpV a )dii 



t JMt 



t JMt 

pV a n Q dp tl , 



(3.21) 



pV a n a dp tl - / pV a n a dp to 

M tl J Mt 

1 /•*!/• „ _ _ _ _ 

T a/3 (0)(V Q V> + V p V a )dp+ I pV a n a dp t0 

M t0 



(3.22) 



or 



E t {4>) + E t (dust) = E ((f>) + E (dust) + 



1 f l 



to JM 



T Q/3 (0)(V Q V^ + V V a )dfi, (3.23) 
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where by definition and Eq. (13.201) . for any t, 

E t (dvsb) = / pV a n a dfi t . (3.24) 
JMt 

We see that the last term in Eq. (13.231) can be zero only when V is a Killing vectorfield. 
We therefore arrive at the following result about the total slice energy with respect to 
the fluid itself. 

Theorem 3.1 The total slice energy with respect to the timelike vectorfield V, tangent 
to the dust timelines, of the scalar field- dust system satisfying the field equations A3.9\) . 
satisfies 

E t {<p + dust) = E {<f) + dust) + ~ f 1 f f a %<p){V a Vp + V p V a )dfi. (3.25) 

Z J t JMt 

In particular the slice energy of the scalar field-dust system is conserved when V is a 
Killing vectorfield of g. 

We also conclude that the property of the conservation of slice energy for dust is a 
conformal invariant. However, when V is not a Killing vectorfield, we see that there 
is a nontrivial contribution to the slice energy coming from the combined effect of the 
stress tensor of the scalar field generated by the conformal transformation coupled to 
the non-stationarity of the spacetime due to the lack of a Killing vector. Note that this 
contribution is also nonzero even in the special case that Eq. (13 . 18[) is assumed for in that 
case the first term in f al3 ((j))(V a V^ + V pV a ) is zero because V a V p {V a Vp + VpV a ) = 0, 
but the whole combination is still not zero as there are additional terms coming from 
the contributions of the other terms in Eq. (12. lip (unless the fluid satisfies an extra 
condition - see below). 

We now proceed to see how this result changes when we assume that (V, g) is filled 
with a perfect fluid. With our conventions the stress tensor of a perfect fluid with energy 
density p and pressure density p is T a/3 = (p + p)V a Vp — pg a /3 and the fluid satisfies the 
field equations 

f'Rap - \g«pf - V Q V^f + g aP U g f = {p + p)V a V p - pg af3 . (3.26) 
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In this case, because the energy-momentum vector P a = pV a , we have P a n a = pV a n a 
and so the slice energy with respect to the timelike vectorfield V in the Einstein frame 



is again 
Then 

Since V a V° 
find 



E t (fluid) 



pV a n a dp t 



(3.27) 



Mt 



V a J fluid — V c 



(p + p)V a Vf 3 + {p + p)V a V a V 13 - V p p. (3.28) 
1, we have V a V pV a = and so on multiplication of Eq. (I3.28P by Vp we 



(p + p)V a 



pV a V a + V a {pV a ). 



(3.29) 



Integrating Eq. (13.291) on the spacetime slab V and using Eq. (13. 2 1 p to re-express the 
p-term in f!3.29j) we have 



-pV a V a dp 



pdp tl 



pdp 



to 



Mt 



(3.30) 



't J Mt Jto J Mt 

and therefore we obtain from (13. 7p the general energy transport equation in the form 

1 ^ 



E t , (d>) - E. 



f af3 {<P){V a Vp + V p V a )dfi 



t JM 



pV a V a dp + R t0 (fluid) - R h (fluid). (3.31) 

'to JMt 

We thus arrive at the following result. 

Theorem 3.2 The total slice energy of the scalar field - perfect fluid system satisfying 
the field equations 113.26}) . depends upon the integrated pressure according to the formula 



^(0+fiuid) = £ to (0+fiuid)+^ / 1 / f aP (<p)iy a Vp+VpV a )dp- [ 1 f pV a V a dp. 

^ Jto JMt Jto JMt 

(3.32) 

In particular, the slice energy is conserved when V is a Killing vectorfield for g (stationary 
spacetime) . 
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When V is not a Killing vectorfield, this slice energy is not generally conserved and this 
is true even in the special case of a fluid with zero expansion, V a V a = 0, for which 
the last term in Eq. f 1 3 . 3 2 j) is zero. In this case the term depending on the scalar field 
continues to have a nonzero contribution to the total slice energy. This term is given by 



T a ?(<p)(y a % + v p v a ) = d a (j ) d^{v a Vp + v p v a ) - v Q v Q (9>9 A - 2V{4>)) (3.33) 

and so Eq. fj3.32j) becomes 

^(0 + fluid) = £ to (0 + fluid) + J r / 9 a 0^0(V a V^ + VpV a )dfl 

1 J to J Mi 

- I 1 ( (ld x <t>d X (/> + p-V(<t>))V a V a dp,. (3.34) 
Jto JMt 1 

We conclude that the only other possible case for which we have slice energy conservation 
is when V is not a Killing vectorfield for g, but the alignment condition (13.181) holds 
(making the middle term in Eq. (13.341) equal to zero) and the fluid has in addition zero 
expansion (last term in Eq. ( 13.341) is zero). 



4 Discussion 

The results of this paper allow us to make some comments concerning the problem 
of deciding which of the two frames (or metrics), Jordan or Einstein, is the physical 
one, meaning in which of the two representations of the dynamics test particles follow 
geodesies (assuming the validity of the principle of equivalence). In the case where test 
particles follow geodesies in both frames, one says that the two conformally related frames 
are physically equivalent. The main results of Section 3, in particular Eq. 13.251 (as well 
as its pressure extension - Eq. I3.32p . were proved under the implicit assumption that the 
vectorfield V a of the dust streamlines generates geodesies in both the Jordan frame and 
the conformally related Einstein frame. 

But does V a always generate a geodesic in the latter frame so that V^V a T^ st = 
V ' a (pV a ) (cf- Eq. 13.161) ? In general, it will not do so and the two frames will not be 
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physically equivalent. In this case, we have 

^V Q T d f st = V a {pV a )%V^ + p(V a V?)V a V^ (4.1) 



so that the result of Theorem 13.11 becomes, 

-ti 



^(0 + dust) = £ o (0 + dust) + ^ I 1 [ f^(<t>)(V a Vp + V p V a )dji 

1 J to JMt 



>t JMt 



t JMt 

h 



E (<f) + dust) + ~ / 1 / (f a %<p) - pV a V p ) (V a V p + VpVaW. 

Z J to JMt V ' 

(4.2) 



We may therefore conclude that the expressions for the total slice energy in the two 
situations considered here, namely, when test particles follow geodesies in both metrics 
g, g (cf. Eq. I3.25j) . or only in the original Jordan frame metric g (cf. Eq. 14.21) . are different 
and in the latter case there is an extra term contributing to the total energy (i.e., the 
dust term in the integrand in the last term in Eq. 14.21) . This additional contribution 
will appear as a measurable quantity which, if measured to be nonzero, will lead us to 
conclude that the two conformally related frames cannot be physically indistinguishable. 
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